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We discuss the finite ma corrections associated with the computation of sea quark matrix elements. We find 
them to differ from the standard normalization used for valence quarks and to depend strongly on the Lorentz 
structure of the current under consideration. Phenomenological implications of these results are briefly discussed 
in two examples. We also mention how the magnitude of the correction factors can be reduced by using a 2-Iink 
improved action. 



Matrix elements of bilinear quark operators, 
measured with the Wilson fermionic action are 
usually remormalized by a factor: 



2k(1 + ma) [1 + 0(ay)] 



(1) 



where, in the context of mean-field improved per- 
turbation theory , 



k — kU{ 







K 



and 



ma = 4 

2k 



(2) 



Here, we would like to show that this formula 
does not apply when the current is part of a "dis- 
connected" quark loop (as is the case for exam- 
ple in computations of sea quark matrix elements 
or in determinations of the mass of flavor singlet 
mesons). Instead, we will find that in these situa- 
tions, the factor {1 + ma) in (1) has to be replaced 
by another ma dependent factor, which depends 
on the Lorentz structure of the current under con- 
sideration. 

The limitations of the (1 -I- ma) factor are eas- 
ily understood by noting the assumptions going 
into its derivation It is obtained by trying 
to match the lattice and continuum propagators 
of a zero-momentum quark on mass shell. One 
can reasonably argue that this is the relevant sit- 
uation for a valence quark inside of a hadron at 
rest (although there is room for improvement at 
small and intermediate values of ma). However, 
this is clearly insufficient to deal with the case of 
closed quark loops where one has an integration 
over 4-momentum. This is even more so given 
that, on the lattice, the integral will be sensitive 



to the contributions coming from the doublers at 
the edges of the Brillouin zone. As is well known, 
Wilson fermions describe 16 species with masses 

ma ma + 2r ma -t- 4r ma + 6r ma + 8r 
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where the degenaracies have been indicated be- 
tween braces. When ma 0, the doublers decou- 
ple (apart from giving rise to the anomaly) and 
continuum results are recovered However, in 
current lattice simulations, the condition ma <^ r 
is not always satisfied, so that the degeneracy 
between fundamental fermions and doublers is 
only partially lifted. This induces lattice artifacts 
which are reminiscent of the naive fermion case: 
closed quark loops involving a scalar current are 
overestimated, whereas those associated with a 
pseudoscalar or an axial current are underesti- 
mated. As we will see below, this requires large 
corrections even for values of ma as low as 0.1 
(at r = 1) which on current lattices corresponds 
roughly to the mass of the strange quark. Dispos- 
ing of a method for estimating the magnitude of 
these renormalizations is therefore of considerable 
phenomenological importance. 

The case of the pseudoscalar current was con- 
sidered several years ago in an interesting series of 
papers by Smit and Vink Q . Their starting point 
is an exact relation (derived from the anomalous 
Ward identity) between the pseudoscalar charge 
(Q5 = j d'^xj^{x)) and the topological charge of 
the background gauge field {Qt). 

mQ5 = Qt (3) 
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By imposing that this continuum relation should 
be satisfied for Wilson fermions on the lattice, 
they were able to deduce the appropriate ma 
dependent renormalization factor for the lattice 
pseudoscalar current. Here we would like to ob- 
tain the corresponding normalization for other 
currents, in cases where an exact relation of the 
type (3) is not available. For this we will make 
the assumption that the ma dependence of the 
renormalization factors is independent of what 
background field is considered, so that it can be 
computed in perturbation theory for a weak ex- 
ternal field. The relevant diagrams at lowest or- 
der are triangle diagrams and by comparing their 
lattice and continuum values, we can deduce the 
appropriate correction factors (see Q for details). 
For example, if we consider the case of the scalar 
current, we would find in the continuum the per- 
turbative relation: 



(4) 



giving the zero momentum matrix element of the 
scalar current in a background gauge field of mo- 
mentum p, where C is a known dimensionlcss 
function. Doing a similar computation on the lat- 
tice, we would find: 



m < ipip >— L{pa, ma)F'^jj[p) 



(5) 



where the lattice dimensionlcss function L can 
now depend separately on pa and ma (but will 
give back C when both tend to 0). For simplic- 
ity, we will consider here only the case of smooth 
background fields (i.e. p — ^ 0). Then, we deduce 
from (4) and (5) that the appropriate ma depen- 
dent correction factor for the scalar current is: 



Ks{ma) = C{0)/L{0,ma) 



(6) 



Similar computations can be done for other cur- 
rents (in particular, for the pseudoscalar current, 
one recovers the result of Smit and Vink men- 
tioned above). The results are given in fig. 1 for 
the scalar (ks), axial {ka) and pseudoscalar (kp) 
currents (to increase the lisibility, the inverse of 
the renormalization factors have been plotted). 
We would like to reiterate that those factors are 
to replace the (1 + ma) part of the renormaliza- 
tion factor whenever "disconnected" quark loops 
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Figure 1. Inverse normalization factors for var- 
ious sea quark currents versus quark mass. The 
valence factor has been included for comparison 
[Ky^ — 1/(1 -I- ma)]. From top to bottom: scalar, 
valence, axial and pseudoscalar 



are involved. They merely correct for the differ- 
ent "response function" of the lattice and contin- 
uum fermionic currents to the same background 
gauge field. They differ from current to current 
because of the different incidence of the doubling 
problem in different channels. The rest of the 
renormalization factor [1 + 0{av)] is as usual. It 
appears that these correction factors are rather 
large and will therefore strongly influence the in- 
terpretation of simulation results. To illustrate 
this, we now discuss their application to two in- 
teresting problems. The first is the computation 
of the mass of the rj' . In a quenched simulation, it 
can be estimated from the ratio of a 2-loop to 1- 
loop amplitude |^ . As a consequence of the above 
discussion, we expect that they will have different 
renormalization factors: The 1-loop part behaves 
as a "valence" quark (and will therefore pick-up 
a (1 -|- ma) correction for each current) whereas 
the 2-loop diagram corresponds to a "sea" quark 
situation (in which case, the pseudoscalar factor 
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of fig. 1 is relevant). In figure 2, we compare the 




Figure 2. rjo mass vs quark mass after renormal- 
ization (squares) and before (plusses). Also in- 
cluded are the chiral extrapolation of the raw data 
(diamond) and the estimate from the Witten- 
Veneziano formula (circle), as obtained in [6]. 



mass of the r]' after correction with the raw lat- 
tice data p. The trends are very different: after 
correction, the mass of the rj' is essentially inde- 
pendent of the mass of the quark up to the strange 
quark mass with a chiral limit that is in much bet- 
ter agreement with the estimate obtained from 
the Witten-Veneziano formula (We assume that 
the drop in the renormalized data at low quark 
masses can be attributed to finite size effects or 
to the zero-mode shift effect Q). The agreement 
between the renormalized data and the Witten- 
Veneziano estimate doesn't come as a surprise, 
since the correction factor that we are advocat- 
ing here is the same as the one used by Smit and 
Vink to compute the topological charge by the 
fermionic method. As a second example, we men- 
tion the contribution of the strange quark to the 
mass of the nucleon. If we replace the (1 + ma) 



correction factor used in Q], with the sea quark 
scalar correction factor from fig. 1, we find that 
the strange quark contribution to the mass of the 
nucleon comes down, from 30% to 21%. 

Since the correction factors given above are 
generally large, there is much interest in try- 
ing to reduce their magnitude by using modi- 
fied actions. We found in ||] that for the axial 
and pseudoscalar current considerable improve- 
ment can be achieved by using the 2-link action 
of Hamber and Wu For the scalar current, 
on the other hand, this kind of action doesn't 
seem to help much and it is probably more in- 
teresting to consider other alternatives, such as 
staggered fermions for example or working with 
naive fermions and dividing the final result by 
16 §. Finally, we would like to recall that the 
renormalization factors presented here have been 
obtained by assuming a lattice simulation dom- 
inated by low momentum gauge fields. In more 
general situations it might be necessary to con- 
sider the "response functions" at finite momen- 
tum This however is left for later studies. 
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